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ABSTRACT 

We argue that a variety of solar data suggest that the activity- cycle timescale variability of the total 
irradiance, is produced by structural adjustments of the solar interior. Assuming these adjustments 
are induced by variations of internal magnetic fields, we use measurements of the total irradiance and 
effective temperature over the period from 1978 to 1992, to infer the magnitude and location of the 
magnetic field. Using an updated stellar evolution model, which includes magnetic fields, we find that 
the observations can be explained by fields whose peak values range from 120k to 2.3k gauss, located 
in the convection zone between O.959i?0 and O.997i?0, respectively. The corresponding maximal radius 
changes, are 17 km when the magnetic field is located at O.959i?0 and 3 km when it is located at 
O.997i?0. At these depths, the W parameter (defined by Alni?/AlnL, where R and L are the radius 
and luminosity) ranges from 0.02 to 0.006. All these predictions are consistent with helioseismology and 
recent measurements carried out by the MDI experiment on SOHO. 

Subject headings: Sun: interior — Sun: magnetic fields 



1. INTRODUCTION 

Direct satellite measurements of total solar irradiance 
(Willson & Hudson 1991; Frohlich & Lean 1998), show 
that it varies almost in phase with the solar activ- 
ity cycle, and that its relative variation in amplitude 
is about 0.1%. Although the most common explana- 
tion of this change, involves the effect of magnetic re- 
gions and network (superposed to an otherwise con- 
stant solar photosphere), an alternative explanation is 
that the cyclic variation is primarily due to structural 
changes of the solar interior. Two arguments support 
this second possibility. Firstly, direct measurements 
of the spectrum of p-mode oscillations, show variations 
(Libbrecht & Woodard 1990) that can be best explained 
in terms of structural changes (i.e. changes in pressure, 
density, radius, etc) (Lydon et al 1996; Antia et al 2000; 
Balmforth et al 1996). Secondly, measurements of the "ef- 
fective temperature" by Gray and Livingston (1997b), also 
vary nearly in phase with the solar cycle. This variation is 
obtained by monitoring the neutral carbon A5380 line in 
the solar flux spectrum. They present several arguments 
supporting their contention that what they measure is the 
gas temperature in the deep photosphere, free from the in- 
fluence of sunspots, faculae and small scale magnetic flux 
tubes. Their main argument rests on their simultaneous 
monitoring of C I A5380, Fe I A5379 and Ti II A5381. 
They use the ratios of spectral line depths, C I A5380 to 
Fe I A5379 and to Ti II A5381, as temperature indicators. 
Since these three lines have different excitation potentials 
(1.57 eV for Ti II A5381, 3.69 eV for Fe I A5379, and 7.68 
eV for C I A5380), the consistency of the results for the 
two line depth ratios, shows that the temperature indi- 
cators arc free from the influence of variations of faculae 
and small scale magnetic flux tubes, which would be ex- 
pected to affect each line in a different way. See Gray & 



Livinston (1997a,b) for the entire discussion. In order for 
the gas temperature to change, the energy flow from the 
interior must also change. Structural adjustment will be 
able to accomplish this. Since the relative radius variation 
A In R is very small (Emilio et al 2000), 

A In L « 4A In T, (1) 

where L and T are luminosity and temperature, respec- 
tively. The temperature change measured by Gray and 
Livingston (1997) is 1.5 ± 0.2 K in one cycle, which cor- 
responds to A In T ~ 0.025% and AlnL ~ 0.1%. This is 
approximately the entire cyclic variation of the total irra- 
diance. 
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The most obvious way the solar structure may change, 
is in response to variations of the internal magnetic field 
during the dynamo build-up and decay. In order to com- 
pute this effect, we introduce magnetic variables as new 
stellar structure variables (Lydon & Sofia 1995), in addi- 
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tion to the conventional ones (Guenther et al 1992). Since 
magnetic field is not a scalar, we have to use at least two 
variables to describe it: the magnetic energy density \ 
(Lydon & Sofia 1995), and the ratio of the magnetic pres- 
sure to the magnetic energy, 7— 1 (Endal et al 1985). This 
method was developed by Lydon and Sofia in 1995, in 
which 7 was treated as a parameter, and it was subse- 
quently applied by Lydon, Guenther and Sofia (1996), to 
explain the observed variation of solar p-mode oscillations 
with the solar cycle (Libbrecht & Woodard 1990). Here 
we want to update this method by: (i) treating both \ and 
7 as new structure variables, as done for \ m Lydon and 
Sofia (1995), (ii) taking into account the influence of mag- 
netic fields on radiative opacities, (iii) taking into account 
all timc-dcpcndcnt contributions to the equations of stellar 
structure (since we want to treat short time scales), (iv) 
modifying the radiative loss assumption of a convective el- 
ement (as discussed in section 3), (v) updating the code 
by using the new version of Yale Stellar Evolution Code 
(YREC7) and (vi) removing the perfect gas approxima- 
tion, which was assumed when calculating some first- and 
second-order derivatives associated with magnetic fields in 
Lydon and Sofia (1995). Because luminosity, radius and 
temperature variations are sensitive to the location, in- 
tensity, orientation, and distribution of the perturbation 
magnetic field B = (Snxp) 1 ^ 2 (p is the gas density), we 
can use the measured yearly-averaged irradiance, temper- 
ature and radius variations, to determine the solar interior 
magnetic field as a function of mass and time. 
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Fig. 2. — The measured solar photospheric temperature varia- 
tions from 1978 to 1992 (Gray & Livingston 1997b) and the yearly 
mean. 

2. SOLAR IRRADIANCE AND EFFECTIVE TEMPERATURE 

RECORDS 

Space-based solar irradiance measurements have been 
reviewed by Frohlich and Lean (1998), and a compos- 
ite solar irradiance record for the period from 1978 
to the present can be downloaded from the website 
http://www.ngdc.noaa.gov/. In our magnetic perturba- 
tion model calculations, we use one year as the timestep, 
since variations on the shorter time scales may be at- 
tributed to the effect of magnetic regions and network. We 
compare the results with the yearly means of total solar 
irradiance, as shown in Figure 1. The standard variances 
range from 0.1 Wm~ 2 to 0.5 Wm~ 2 . 

The spectroscopic temperature variations of the sun, 



measured by Gray and Livingston (1997) over the period 
from 1978 to 1992, are shown in Figure 2 (reproduced from 
their Figure 9). The points are the yearly means, the zero 
point being chosen arbitrarily. 

3. METHOD 

3.1. Definition of variables 

Standard solar models (SSM) use pressure, temperature, 
radius and luminosity as the structure variables, and mass 
M r interior to a radius r, as the independent variable. 
In order to self-consistently take into account magnetic 
effects, Lydon and Sofia (1995) introduced two new mag- 
netic structure variables to model the sun: the magnetic 
energy per unit mass \-> an d the effective ratio of spe- 
cific heats for the magnetic perturbation, 7. The former 
describes the magnetic perturbation strength, and the lat- 
ter describes the tensor feature of the magnetic pressure. 
In general, the determination of x an d 7 requires a com- 
prehensive understanding of turbulent dynamics in the so- 
lar convective zone (Brummell et al 1995), an undertaking 
that is impractical at present. Therefore we take x an d 7 
as free parameters (functions) that need to be determined 
at this stage, instead of using dynamo equations to deter- 
mine them. 

Our aim is to use the measured solar irradiance and ef- 
fective temperature variations given in Figures 1 and 2 to 
determine x(^D,t) and ^(Mmt), where Mjj and t are 
mass depth and time, respectively. Using x(Mr>,t), we 
can obtain the solar interior magnetic field B(Mrj,t) = 
(87TX/3) 1 / 2 . Since these observations are not enough to de- 
termine the two new variables uniquely, we assume they 
have a gaussian-like shape, and also utilize helioseismic 
((Antia et al 2000)) and radius ((Emilio et al 2000)) ob- 
servations, as additional constraints. The mass depth Mrj 
is defined as 

Md = log 10 (l — M r /M ). (2) 
This replaces the mass variable M r when we want to de- 
scribe the location, orientation and form of all magnetic 
perturbations. 

For our purpose we can use the toroidal (B t ) and 
poloidal (Bp) components to express a magnetic field vec- 
tor B = (B t ,B p ). The magnetic energy density variable x 
can be written as, 

X = (B 2 /8tt)/p, (3) 
where B = (£> 2 +B 2 ) 1 / 2 is the magnitude of the magnetic 
field vector. The magnetic field direction variable 7 can 
be written as, 

1 = (2B 2 t +B 2 p )/B 2 . (4) 

The magnetic pressure P x can be defined in terms of the 
new magnetic variables as 

P x = h- 1)XP- (5) 
The value of 7 depends on the field geometry. Parallel to 
the field lines, P x is zero, so 7 = 1, whereas perpendicular 
to the field lines, P x is a maximum, so 7 = 2. 

3.2. Consequences on the equations of stellar structure 

The continuity equation remains the same. The hy- 
drostatic equation is replaced by the equation of motion, 
which including the Lorentz force, is 

d 2 r GM r „ 1 , s . -. 

= —Pr + ^(V x B) x B, 
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where G is the gravitational constant. The radial compo- 
nent is 



d 2 r 



dP GM r 

-ri> 



<9P V 



dr r 2 dr 

The magnetic field has both a magnetic tension and a 
magnetic pressure. This can be seen by writing the 
Lorentz force (l/47r)(V x B) x B, as the sum of the 
magnetic tension (B ■ V)P/47r and the magnetic pres- 
sure — V(B 2 /8ir). The former is anisotropic, and the 
latter is isotropic. Both contribute to the last term of 
the radial motion equation, and are incorporated into P x . 
For example, for a purely radial magnetic field, the mag- 
netic tension force, {1/ Air)B r dB r / dr = (l/%n)dB 2 /dr = 
(l/8ir)dB 2 /dr, which cancels the magnetic pressure force. 
As a result, P x = in this case. 
Defining 

Pt = P + Pv 



as the total pressure, the equation of radial motion be- 
comes 

dP T GM r 1 d 2 r 



dM r 4-irr 4 Airr 2 dt 2 ' 
where the last term represents the inertial force. Since 

P = p{Pr,T,x,i), 
the equation of state becomes 



(6) 



where 



dp 
P 
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dP T 
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dlnp 
d\nP T 
d\np 

~dlnP 
<91np 

dhix 
dlnp 
d In 7 



dT d X ,d 1 
d— — v v — , 

T X 7 



at constant T, x, 7 
at constant Pt, X-> 7 
at constant Pt, T, 7 
at constant Pt, T, x 



Since we use the one-dimensional stellar evolution model, 
we can not model the transverse component of the equa- 
tion of motion. In order to handle this transverse compo- 
nent, we need at least a two-dimensional stellar evolution 
model. Since such a model is not available yet, we only 
consider the radial motion equation in this paper. 
The energy conservation equation becomes 



dL 



,dS 7 



1 du 
p~~di' 



(7) 



dM r dt 
Here 

TdS T = dQ T = dU + PdV + d x 

= dU T + P T dV - (7 - l)(x/V)dV 

,rr, s ( P T Sv\ , Pt5v' , 

= CpdT - -dP T + 1 + — — d X + — di%) 

P \ pax J P°H 

is the first law of thermodynamics including magnetic 
fields. The total internal energy Ut = U + x and the total 
entropy St = S + x/P. c p is the specific heat at constant 
pressure. If 7 is constant, then dj = 0, and equation (8) 
reduces to equation (75) in Lydon & Sofia (1995). The 
symbol u — aT 4 , is the radiation energy density, where a 



is the radiation constant, du/dt appears in the full energy 
conservation equation, 

du „ ~ , ^dSr , 



lH +V - F = ^- T -dt 



F = uv is the radiation energy flux vector, and v is the 
photon diffusion velocity. 

The equation of transport of energy by radiation, 

^L = ? Kj_ 

dM r Mir 2 acr 4 T3' V ' 

does not change in form, but the magnetic field affects the 
opacity coefficient n. Here c is the speed of light. See 
appendix A for an estimate of this effect. 

The equation of energy transport by convection, 
dT PGM. 

MP=-P4^ V < (10) 
does not change form either, but the convection temper- 
ature gradient, V, with a magnetic field, is different from 
that without it. Lydon and Sofia (1995) have demon- 
strated how to account for magnetic effects in the mixing 
length theory, and hence obtain V, provided 7 is constant. 
Their method can also be used to derive V when 7 is vari- 
able, as we are about to do in this paper. See appendix B 
for the concrete expression for V in our case. 

3.3. Numerical implementation 

The stellar evolution code solves the stellar structure 
equations (such as Eqs. [6], [7] and [9] or [10], and the 
continuity equation, dr/dM r — l/47rr 2 p), with suitable 
initial and boundary conditions. In the model, M r and t 
are the independent variables, and the conventional struc- 
ture variables are Pt, T, r and L, while the magnetic 
variables are x an d 7- 

To follow the yearly variations of solar global parame- 
ters with a stellar evolution code, requires a timestep no 
larger than one year. Such a small timestep requires a very 
precise code. For example, because Lq has increased by 
about 30% during the lifetime of the sun, the relative mean 
rate of change is about 10~ 10 per year. To achieve such 
a high precision is a challenge. Fortunately, Yale stellar 
evolution code (YREC) meets this need, as shown before, 
(Lydon & Sofia 1995). We follow Lydon and Sofia (1995) 
by modifying YREC7, a new version released in May 1999, 
in order to accommodate the magnetic effects described 
above. The reason why YREC permits a small timestep 
and achieves such a high numerical precision, is because it 
uses analytical formulae (Prather 1976), rather than nu- 
merical methods, to solve the linearized stellar structure 
equations. 

In order to accommodate the various magnetic effects 
described above, we must first write a routine to specify 
X and 7. We assume x — Xm(t)F{M]j). The maximum 
magnetic energy density Xmit), is to be determined from 
solar activity indices. The yearly-averaged sunspot num- 
ber, Rz, is the most widely used solar activity index. From 
numerical experiments, we find that the results are sensi- 
tive to the function form of B m on R z . If the maximum 
magnetic field in the solar interior, B m , is related to Rz 
via 

B m = P o {190 + [1 + log 10 (l + Rz)] 5 }, (11) 

then by adjusting B , we can nearly match the measured 
cyclic variations of irradiance and effective temperature. 
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The reason why such a functional form of B m is cho- 
sen is that B ~ 20 kG, at a depth of M D = -4.25 
(r = O.96i? ) in 1996. This result is inferred from he- 
lioscismology (Antia et al 2000) when Rz was at a mini- 
mum. Using this prescription for B mi the value of B m is 
about twice as large at the maximum of solar cycle, as it 
is at the minimum. 

F(Mjj) specifies the distribution of \, an d is to be de- 
termined by fitting the measured irradiance and effective 
temperature variations. F(Mjj) has infinite degrees of 
freedom and thus cannot be determined uniquely by ob- 
servational results, which have finite degrees of freedom. 
However, we can remove this degeneracy by assuming a 
field shape of the form, 

F(M D ; M Dc , a) = cxp[-i(M D - M Dc ) 2 /a 2 ], (12) 

where Md c specifies the location and a specifies its width. 
This gaussian profile allows us to pinpoint the location of 
the required magnetic field, by using observations of cyclic 
variations of irradiance and effective temperature. 
We either fix 7, or use 



7 (M D ) = 1 + (B/ Bm) 1 / 5 (R z / 200), 



(13) 



to express the spatial and temporal variations of the direc- 
tion of magnetic fields. This implies that the more intense 
the magnetic field, the larger its toroidal component. The 
profile of 7 is arbitrary, but its value must lie between 1 
and 2. Although the detailed results are sensitive to the 
functional form of 7, the qualitative features are same. We 
find that the magnetic effects are maximized when 7 = 2, 
so we set 7 = 2 in subsequent computations. The main 
effect of a smaller 7, is to increase the magnetic field re- 
quired to reproduce the observed variation of the sun. 

As the magnetic fields that cause the solar activity cy- 
cle are believed to exist outside of the solar core (where 
the thermonuclear reactions occur), we do not consider the 
influence of magnetic fields on the energy generation rate. 

Magnetic fields affect the equation of state in the fol- 
lowing way. If (3, the ratio of the gas pressure to the total 
pressure when magnetic fields are present, is defined as 



1 - 



\«T* (7"1)XP 



(14) 



P T Pt 

(a is the radiation constant), then the density of a gas can 
be determined from 



f3P T = 1lpT(l + E)/p a . 



(15) 



TZ is the gas constant, E is the number of free electrons per 
nucleus (determined by solving the Saha equation) and p, a 
is the mean atomic weight per atom (Prather 1976). Note, 
E also depends on (3. In the stardard case, since x = 0, /? 
does not depend on p. But in the nonstandard case dis- 
cussed in this paper, (3 depends on p. Consequently, not 
only p, but also its first derivatives a, 5, v, and v' , and its 
second derivatives with respect to T and Pt, all need to 
be calculated by iteration (which is tedious). In the orig- 
inal implementation of this method by Lydon and Sofia 
(1995). The second derivatives were approximated by as- 
suming a perfect gas, which needs no iteration. In this 
code upgrade, we modify the equation of state routine so 
that we can accurately calculate all these thermodynamic 
variables by iteration. Although this improves the numer- 
ical precision, it does not make any qualitative difference. 



V x and V 7 (defined in Appendix B) are also calculated by 
numerical derivatives, as are those time derivatives which 
appear in the equations of stellar structure. 

The qualitative difference comes from /' (see Appendix 
B), which represents the influence of magnetic field on ra- 
diative loss of a convection element. The convection tem- 
perature gradient is now given by 

V = V ad + (y'/V'^C)(l + y'/V) - A m , (16) 

where 

A m = f'[{v/a)V x + (*//a)V A ]V ad . (17) 

All symbols used here are defined in Appendix B. When 
/' = 1 and 7 is fixed, equation (17) reduces to 

A m = (^/a)V x V ad . (18) 

However, under those assumptions, the calculated 
cyclic variation of effective temperature is in antiphase 
with the solar activity, contrary to the observations 
(Gray & Livingston 1997b). In order to get a proper fit 
to the observations, it is necessary to assume /' = 3, a 
value that will be used in all our subsequent calculations. 

We use Eq. (A2) to compute the global magnetic effect 
on the radiative opacity, and find that this has little influ- 
ence on the effective temperature. On the other hand, the 
need of /' = 3 implies that the energy loss of a convective 
element due to turbulence associated with the local mag- 
netic field, is much more efficient than the radiative loss 
enhancement of the convective element due to the global 
magnetic field. 
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Fig. 3. — The difference between the temperature gradient and 
the adiabatic gradient as a function of the mass depth in the so- 
lar interior. The dashed line marks the location of lower boundary 
surface of the convection zone. 

We do not redistribute the mass grid points, nor do we 
predict stellar structure variables for the next time step, 
when the timestep is smaller than a solar cycle period, 
since these operations will generate an error much larger 
than 10~ 10 per year for luminosity. All convergence cri- 
teria for iteration are tightened by requiring the relative 
error to be no larger than 6 x 10~ 16 or so. In this modified 
code, the matching point is at Mo ~ —8.6. 

The conventional model parameters are, the mixing 
length parameter (the ratio of the mixing length to the 
pressure scale height, usually denoted as a parameter), the 
initial hellium mass abundance (Yq), and the initial heavy 
element mass abundance (Z). The magnetic parameters 
introduced in this paper are the peak mass depth Md c 
of the magnetic energy density \, the width of \ (^■ e -,< J ), 
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the magnetic strength parameter Bq, and the turbulent 
loss parameter of energy of a convective element /'. It 
is well known that the conventional parameters can be 
determined by present observations of solar radius and lu- 
minosity. Through numerical experiments (i.e., model cal- 
culations), we find that the magnetic parameters (Md c , 
<t, _B and /'), can be determined by the observed cyclic 
variation of solar radius, luminosity and effective temper- 
ature. We fix /' = 3 in all subsequent runs, as required by 
the observed variation of effective temperature. 7 = 2 is 
fixed to maximize the magnetic effects. In order to probe 
the observation-required location, strength and profile of 
the solar interior magnetic field, we change Md c , o and 
Bq from run to run. 

4. CONSTRAINTS ON THE SOLAR INTERIOR MAGNETIC 

FIELD 

4.1. The overshoot layer 

At present, it is generally believed that the magnetic 
field is stored only in the subadiabatically stratified over- 
shoot layer in the form of flux tubes, which may break 
loose due to kink instabilities, when a threshold of 10 5 
gauss is exceeded (Caligar et al 1995; Caligar et al 1998). 
Figure 3 depicts V — V a( j vs Mb, in which the subadi- 
abatically stratified overshoot layer approximately corre- 
sponds to the mass depth range: —1.68 < Md < —1.2. 
Mdc = —1.68 marks the base of the convection zone of 
our model sun. Therefore, the maximum yearly-averaged 
field strength B m in the overshoot layer during the so- 
lar cycle should be smaller than the threshold. The field 
with Mdc = —1.45 and a < 0.2, can be considered to 
be confined in the overshoot layer. However, even if the 
threshold value of B is used, the resultant irradiance and 
effective temperature variations are negligible. 

Only when we abandon this threshold constraint (e.g., 
using a field larger than 2.6 x 10 6 gauss), can we get irradi- 
ance variations comparable to the observations, in magni- 
tude and phase; however, even then, the effective temper- 
ature variation is much smaller and in antiphase with the 
solar cycle. Besides, the resulting relative radius variation 
is larger than 0.03%. Since the observations reveal that 
both irradiance and effective temperature variations are in 
phase with the solar cycle, and the relative radius variation 
should be much smaller than 0.01% (Emilio et al 2000), we 
conclude that the magnetic field that gives rise to the ob- 
served variations of irradiance and effective temperature 
(shown in Figures 1 and 2) cannot be completely confined 
in the overshoot layer. 

4.2. The convection zone 

The above arguments lead us to investigate the possi- 
bility that the magnetic field is located totally within the 
convection zone. We set a — 0.5 in this region. In order to 
reproduce the observed cyclic variation of irradiance and 
effective temperature, B will increase if a decreases, and 
vice versa. We choose Md c from the interval (-8.6,-1.68], 
which covers the convection zone of our numerical solar 
model. For each selected value for Md c , we change Bo 
in order to reproduce the observed cyclic variations of so- 
lar irradiance and effective temperature by using our solar 
evolution code. 



4.2.1. The lower part 

We find that the lower part of the convection zone, 
— 1.68 > Mdc > —4.2, is ruled out because it produces 
radius changes that are too large (Ami? > 0.01%) and/or 
because it requires a too strong magnetic field (B > 5 x 10 4 
G). 

4.2.2. The intermediate part 

For the next region up, —4.2 > Md c > —7.8, or 
0.959i? Q < R < 0.997i? Q , we find that our model can re- 
produce the observed cyclic variations of irradiance and ef- 
fective temperature by using a possible equipartition mag- 
netic field, B m < 5 x xlO 4 G. In this case, the am- 
plitude of the radius changes is from 17 to 3 km, and 
W = A In R/ A In L ranges from 0.02 to 0.006. 

In order to pinpoint the location of the required mag- 
netic fields in the allowed region of solar convection zone 
(Mdc £ [—7.8, —4.2]), we need more information. We ob- 
tain the allowed region by using the observational infor- 
mation for solar irradiance and effective temperature cyclic 
variations, and by assuming the corresponding solar radius 
variation A In R to be much smaller than 0.01%. If we have 
the direct observational information for A In R or the W 
parameter, we can pinpoint the location of solar interior 
magnetic fields. Since we do not have simultaneous obser- 
vational information for the radius cyclic variation in the 
period ranging from 1978 to 1992, we must invoke some 
physical considerations. 



12 






10 






8 






X 






1 6 
i> 




A ; 


1 4 







2 














-2 -4 -6 -8 

Mass Depth 

Fig. 4. — The ratio of radiative to adiabatic temperature gradi- 
ents as a function of mass depth in the solar interior. The dashed 
line marks the location of lower boundary surface of the convection 
zone. 

Figure 4 shows the ratio of radiative to adiabatic tem- 
perature gradients as a function of mass depth in the sun, 
from which we can see that the layer with Md ranging 
from —4 to —8 is the most unstable convective region. 
In this layer superstrong plasma turbulence is inevitable, 
since the radiative gradient V ra d is much larger than the 
adiabatic temperature gradient V a d- Therefore, this may 
suggest a turbulence generation mechanism for solar mag- 
netic fields. Strong plasma turbulence can generate strong 
small-scale magnetic fields, which form small magnetic 
cells with almost random orientation. The residual field of 
the previous cycle and differential rotation, tend to align 
and bundle these random magnetic "needles" to form flux 
tubes. 
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The solid curve in Figure 5 (choosing Md c = —6.5, 
a = 0.5, -Bo = 7 G) shows the corresponding magnetic 
field distribution in the solar interior in 1989, near solar 
maximum. Figure 6 compares the calculated (solid curve) 
and measured (dash-dotted curve) cyclic irradiance varia- 
tion, while Figure 7 compares the calculated (solid curve) 
and measured (dash-dotted curve) cyclic effective temper- 
ature variation. The irradiance fit (x 2 = 2) is better than 
the effective temperature fit (x 2 = 47). The computed 
radius variation is 5 x 10~ 6 , and the W parameter is al- 
most constant (about 4 x 10~ 3 ). Figure 8 shows the calcu- 
lated structure change generated by the field in 1989 (solid 
curve) . 
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Fig. 5. — Three possible distributions of inferred magnetic field 
in the solar interior in 1989 according to the measured irradiance 
and photosphcric temperature cyclic variations given in Figures 1 
and 2. The vertical line indicates the base of the convection zone. 

Although we can relate the magnetic field gener- 
ation to the plasma turbulence in this region, the 
magnetic field maintenance here faces the old prob- 
lem of magnetic buoyancy. How can we get around 
this problem? Observations and simulations sug- 
gest that the actual dynamics within solar convec- 
tion zone is extremely intricate (Brummcll ct al 1995; 
Kim & Chan 1998; Nordlund 1999). The velocities and 
magnetic fields are complex, exhibiting large-scale struc- 
ture and ordered behavior amidst rapidly varying and 
intense small-scale turbulence. In fact, the numer- 
ical simulations of the solar outer convection zone 
(Kim & Chan 1998; Nordlund 1999) indicate a major 
presence of downward moving plumes. These ordered 
downdrafts may gather small magnetic flux tubes, gen- 
erated in the extremely unstable layer, to form larger flux 
tubes and carry them to the deeper layer. These down- 
drafts may also push down the magnetic flux tubes to bal- 
ance the magnetic buoyancy to form a magnetized layer 
below the most unstable convective region. 

One possibility is that the downdrafts gather and carry 
magnetic flux tubes, generated by turbulence, to depths 
in the convective envelope until some sort of equipartition 
is reached. In fact, Antia, Chitre, and Thompson (2000) 
have explored this possibility. They employ the observed 
splittings of solar oscillation frequencies to separate the ef- 
fects of interior solar rotation, and to estimate the contri- 
bution from a large-scale magnetic field. After subtracting 
out the estimated contribution from rotation, there is some 
residual signal in the even splitting coefficients. This may 



be explained by a magnetic field of approximately 20 kG 
strength located at a depth of M D = -4.25 (r = 0.96i? Q ) 
in 1996. Since the density near Mq = —4.25, is of order 
4 x 10~ 3 , and the downward velocity for the plumes is of 
order 5 x 10 4 cm s _1 , the estimated dynamical pressure of 
the plumes, pv 2 , is equal to or larger then 10 7 dyne cm~ 2 . 
The size of pv 2 is comparable with the magnetic pressure, 
B 2 /8tt, corresponding to a field strength of 20 — 30 kG. 
This demonstrates that a stable magnetized layer in the 
convection zone proper, may form when the complexity 
of convection motion is taken into account. If we wish to 
reproduce the observed temperature and irradiance varia- 
tions by the magnetic field at the depth indicated by he- 
lioscismology, we find that B m ranges from 20 kG to 47 
kG during a solar cycle (choosing _B = 90 G). There- 
fore, in Figures 5-8 we also show this case (dotted curves: 
M Dc = -4.25, ct = 0.5). The irradiance fit (x 2 = 1) is also 
better than the effective temperature fit (x 2 = 43), the W 
parameter is equal to about 2 x 10~ 2 , and the predicted 
cyclic variation of solar radius is equal to about 2 x 10~ 5 , 
which can be tested by measuring the W parameter. 

In fact, using MDI/SOHO data obtained between April 
19, 1996 and June 24, 1998, Emilio et al (2000) found 
W < 2 x 10~ 2 , which is consistent with the above predic- 
tion. This shows that the second case is in agreement with 
all relevant precise observations, including solar irradiance, 
effective temperature, radius, and p-mode oscillation ob- 
servations. The first case is ruled out by the W parameter 
inferred from the observation. 

4.2.3. The upper part 

The upper part of the convection zone, Md c < —7.8, 
is ruled out because the resulting irradiance and effective 
temperature variations are in antiphase with the solar cy- 
cle, which is in conflict with the observations. 

In all cases, to produce the observed luminosity and ef- 
fective temperature variations, the required magnetic field 
must increase with depth. The resulting W parameter also 
increases with depth. This is in agreement with the early 
studies (Dappen 1983; Endal et al 1985). 

4.3. The extended layer 

Another possibility is that these magnetic flux tubes are 
carried downwards into the stable subadiabatically strati- 
fied overshoot layer, to provide a seed field for the dynamo 
operating at the base of the convection zone. Therefore, 
the magnetic field extends from the convection zone to the 
overshoot layer. The following magnetic field distribution 
mimics this case, 

f M700) M D > M Dc 

F = < h(l) - 0.325(M D - M Dc ) M d < M Dc 

[ h(l) + 1.3exp[-2(M D + 5.55) 2 ] M D < -5.55 

(19) 

where h(m) = (M D /M Dc ) m cxp[~m(M D /M Dc - 1)]. 
Choosing Md c = —1.55 in equation (19) and Bq — 250 
G in equation (11), we also get a good fit to Figures 1 
(x 2 — 1) an d 2 (x 2 = 40). The magnetic field profile is 
depicted as a dashed curve in Figure 5. This case is also 
represented by dashed curves in Figures 6-8. The pre- 
dicted W parameter is about 2 x 10~ 2 . It should be noted 
that what contributes to the observed cyclic variations of 
solar irradiance and effective temperature in this case, is 
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the field located within the convection zone proper, while 
the contribution of the field confined in the overshoot layer 
is negligible. MD1/SOHO data do not rule out this possi- 
bility. 
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Fig. 6. — Comparison between the measured (dot-dashed curve) 
and calculated solar irradiance variations. The line style for the 
calculated is the same as that in Fig. 5. 
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Fig. 7. — Comparison between the measured (dot-dashed curve) 
and calculated solar photospheric temperature variations. The line 
style for the calculated is the same as that in Fig. 5. 

5. CONCLUSIONS 

It is possible to locate solar interior magnetic fields using 
the observed cyclic variations of three global solar param- 
eters such as luminosity, temperature and radius at the 
surface of the sun. This provides an alternative to helio- 
seismology as a probe of the solar interior magnetic fields. 



Simultaneous measurements of solar total irradiance and 
effective temperature, can only select an allowed range of 
solar internal magnetic fields, which in terms of magnitude 
and location are consistent with helioseismic observations 
and recent MDI experiment on SOHO. 

Although the observed cyclic variations of solar irradi- 
ance, effective temperature, radius and p-mode oscillation 
frequencies, require a magnetic field component between 
20 and 47 kG, peaked at r = O.96i?0 (within the convec- 
tion zone proper), a stronger component of about 300 kG 
buried in the overshoot layer beneath the base of the con- 
vection zone cannot be ruled out, since the contribution of 
the latter to those observations is negligible. 
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Fig. 8. — The structural changes caused by the magnetic field 
distributions given in Figure 5: relative pressure, temperature, ra- 
dius and luminosity changes from top to bottom. The vertical line 
indicates the base of the convection zone. The line style for the 
calculated is the same as that in Fig. 5. 
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APPENDIX 

AN ESTIMATE OF MAGNETIC EFFECTS ON THE RADIATIVE OPACITY 

A magnetic field affects the absorption processes associated with free electrons, such as electron scattering absorption 
K sc and free-free transition absorption ks- We can estimate this effect in terms of the mean free paths with and without 
a magnetic field, Ib and l 

Kco/KcB cx (l B /l ) 2 w 1 + i( 7 - l)T e 2 f>2 (Al) 

Here 1 /r e is the collision frequency between electrons, n e = k sc + Kg is the absorption component that will be affected 
by magnetic field and f2 e is the electron cyclotron frequency. Magnetic fields also affect the bound-bound absorption, but 
this effect cannot be estimated easily by means of the classic method. Nevertheless, the quantum effect must be much 
smaller than the classic one described above. 
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The influence of magnetic fields on the radiative opacity can only be treated approximately, since we only use opacity 
tables. We assume that we can decompose the total opacity coefficient k into two parts: k = K e + Ki, since n sc and Kg 
do not or only weakly depend on frequency (before taking the Rosseland mean). We use the opacity subroutine provided 
by YREC7 to calculate k, and use equations (17.2) and (17.5) in Kippcnhahn and Weigert (1990) to calculate K e . The 
opacity corrected by magnetic fields, n' , can be expressed approximately as follows 

(7 - 1)t 2 Q 2 

K ' = K_Ke 4TFip!' (A2) 

The approximation originates from the fact that the absorption due to free-free transitions depends on frequency as 

CONVECTION TEMPERATURE GRADIENT 

When both \ an d 7 are considered to be variables, the convection temperature gradient V = (d In T/8Pt) s (s stands 
for surroundings) can be expressed as follows 

V = V ad + (y'/V'^C)(l + y'/V) - A m , (Bl) 

where the initial \ of the convection element is assumed to remain frozen in the surrounding material. V ad = (d In T/d In Pf) 5 
(S stands for entropy) is the adiabatic gradient and y' is obtained by solving the cubic algebraic equation, 

70, C, V, A m and A arc defined by 



= 2Ay' 3 + V'y' 2 + V' 2 y' - V . 


(B2) 


70 - [( C pp)/(2a C T 3 )][(l + (l/3)w 2 M 


(B3) 


C = (g/l 2 m 5)/8H p , 


(B4) 


A m = f'Kv/a)V x + (i//o)V A ]Vad, 


(B5) 


y' = i/[ 7 oC 1 / 2 (v rad -v ad + yi m ) 1 / 2 ], 


(B6) 


A=(9/8)[lo 2 /(3 + cu 2 )}. 


(B7) 



10 = Kpl m , g is the gravity acceleration, H p is the pressure scale height, l m is the mixing length and V rad is the radiative 
temperature gradient. /' is a dimensionless parameter that determines the influence of magnetic field on radiative loss 
of a convective element, V x = (<91nx/<91nPT) s , and V 7 = (dlwy/dlnP^s- In general, magnetic fields tend to inhibit 
convection. 
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